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Abstract 

We perform a symmetry analysis of modified Brans-Dicke cosmo- 
logical equations and present exact solutions. We discuss how the 
solutions may help to build models of cosmology where, for the early 
universe, the expansion is linear and the equation of state just changes 
the expansion velocity but not the linearity. For the late universe the 
expansion is exponential and the effect of the equation of state on the 
rate of expansion is just to change the constant Hubble parameter. 

1 Introduction 

The standard model of cosmology [1] has undergone several modifications in 
the past twenty years. With the discovery [2] of the cosmic microwave back- 
ground radiation the early cosmological models [3, 4] developed into the first 
standard model which had the early radiation dominated and the late dust 
dominated stages. Observational and theoretical consistency has forced mod- 
ifications on the standard model. We now believe that the universe is born 
in a radiation dominated stage, then undergoes an inflationary exponentially 
expanding stage, then becomes radiation dominated again, and then contin- 
ues with matter domination, which by today has evolved into a dark energy 
dominated, slowly but exponentially expanding, stage. The main reason for 
such a complicated history of domination is that according to Einstein cos- 
mological field equations the rate of expansion of the universe depends on 
the equation of state of the matter-energy that fills it. One immediate ques- 
tion which arises is whether there is any consistent modification of Einstein's 
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equations such that the expansion of the universe is independent of its con- 
tent. In this paper we would like to study such a model. The model consists 
of a modified Brans-Dicke-Jordan-Thirry [5] model where the signs of the 
kinetic and potential terms of the "scalar field" are negative when written to 
the right hand (wooden) side of the Einstein equations but is positive just as 
the geometric term containing the square of the time derivative of the scale 
size of the universe when written to the left hand (marble side) . Thus the 
Friedman equation 

is changed into 

o fa 2 k\ n 2 Au 

where is a new geometric field inspired by the Brans-Dicke-Jordan-Thirry 
model. The reason we call <fi a geometric field is that its appearance is quite 
similar to the scale size a of the universe. The terms . . . will be determined 
by the covariant action. In fact, the analogous equation for a Kaluza-Klein 
theory with n internal dimensions is [6] 

3 @ + ?) + r ( "-4 + - = 8 ^ < L3 > 

where b is the scale size of the internal manifold. 

In section 2 we present the basic equations of the model as applied to 
cosmology. In section 3 we perform a symmetry analysis based on dilatational 
symmetries and show that symmetries exist under the conditions that the 
mass of the "scalar field" vanishes and/or the curvature of spacelike sections 
is zero. The solutions are analyzed in sections 4, 5 and 6. 

In section 7 we conclude with a model of expansion for the universe where, 
in the early stages, the mass term of the "scalar field" can be neglected (m = 
0) but the curvature term is important. For closed (k = 1) spacelike sections 
the universe expands linearly. This linear expansion is independent of the 
equation of state. However, the rate of change of the newtonian gravitational 
constant depends on the equation of state. This is in a strong contrast to 
standard general relativistic cosmology where the rate of expansion strongly 
depends on the equation of state. 

Whereas the curvature term is important in the early stages, we argue 
that in later stages, after the universe has expanded to a large size, it can be 
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neglected. The mass term then causes a slow exponential expansion which 
can be identified with the dark energy phenomenon of the present day uni- 
verse. This phenomenon is also independent of the equation of state which 
only slightly influences the exponential rate of expansion. 



2 Basic equations of the model 

The action is the following: 



(2.1) 



where represents the Brans-Dicke scalar field, oo denotes the dimensionless 
Brans-Dicke parameter taken to be much larger than 1, to >> 1. The scalar 
field has a kinetic term —^^d^d^ with the wrong sign and the potential 
of the scalar field contains only a mass term with the wrong sign |m 2 2 . Lm, 
on the other hand, is the matter Lagrangian such that the scalar field does 
not couple with it. R is the Ricci scalar. For simplicity we also restrict our 
analysis to the Robertson Walker metric to emphasize that is necessarily 
spatially homogeneous: 



ds 2 = dt 2 - a 2 (t) 



dx 
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where k is the curvature parameter with k — — 1, 0, 1 corresponding to open, 
flat, closed universes respectively and a (t) is the scale factor of the universe. 
After applying the variational procedure to the action and assuming = (t) 
and energy momentum tensor of matter and radiation excluding is in the 
perfect fluid form of Ttf = diag (p, — p, —p, —p) where p is the energy density 
and p is the pressure and also noting that the right hand side of the 
equation must be zero in accordance with our previous argument on L M 
being independent of 0, the field equations reduce to (dots denote d/dt) 



a , 
+ 3-0 + 
a 
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m + ^r\- + — + — 



= 0. (2.5) 



The following identity shows that the continuity equation could be considered 
as a consequence of these three equations: 



-^(2.3) + 3- x ((2.3) + (2.4)) - x (2.5) 



= p + 3-{p + p) = 0. (2.6) 

(X 



Instead, we prefer to consider more complicated pressure equation (2.4) as 

an algebraic consequence of i (2-5), (2.3) and the continuity equation 

in (2.6). 

We assume the power law for the density 

p = Ca u p = vpH (2.7) 

which implies that the continuity equation (2.6) becomes the equation of 
state: 

so that for dust we have v = —3, p = and for radiation v = —4, p — p/3. 
Now define new unknowns F(a) and H(a) by the formulas 

F=^, H = - - = aHH' + H 2 , — = aHF' + F 2 , (2.9) 

a a 

where prime denotes derivative of a function of a and the relations H = 
H'(a)a = aHH', F = F'(a)a = aHF' were used. Equations (2.3)-(2.5) 
become the equations for unknowns H(a) and F(a) 

tti ^ttt, 2cj , 2 2\ k /4aA p /4uA Ca v /n 

* + + T (F + m ) + * = (t J ? s (t J ^ (2 ' 10) 

- a// (if' + F') + H 2 + - tfF + - (2 - cu)F 2 + — + — m 2 
3 V ! 3 3 V ; 3a 2 3 

aH (\h'+?- F') + H 2 + coHF+^ (F 2 + m 2 ) + = 0. (2.12) 
\2 6 J 6 la A 
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The more lengthy equation (2.11) is still an algebraic consequence of > 
(2.12), (2.10) and the continuity equation in (2.6) due to the linear relation: 

d ^ 2 " 1 + 2F x (2.10) + 3H x ((2.10) - (2.11)) +4Fx (2.12) = 0, 
dt 

so we choose (2.10) and (2.12) together with the state equation (2.8) as 
independent dynamic equations. However, these equations do not form a 
closed system for determining H(a) and F(a) because of <fi in the right-hand 
side of (2.10). To eliminate from the p-equation (2.10), we differentiate 
(2.10) with respect to time, use 

a = aH, <j) = F(j), H = aHH'(a), F = aHF'(a) 

and eliminate (f> with the aid of (2.10) with the result 

(H 2 + HF) H' + (h 2 + - uHF J F' = — H 3 + ^ - ^ H 2 F (2.13) 
V 3 J (x 



+ (^-6) HF2 _2u p , + k 
3a 3a a 6 



v^ + 2) 



H — F 



,2, 



Solving algebraically equations (2.13) and (2.12) with respect to H' and F', 
we obtain the closed system of two first order equations in normal form that 
determines two unknown functions H(a) and F(a): 

(2iu-3)aH ( ^- = (vlu + Q)H 2 + 2(is + 4)uHF+ — [(u + Q)lu - 3]F 2 
da 3 



+ — (uco + 3)m 2 + — [(i/ + 2V + 3], (2.14) 

o (X 

-(2to-3)aH— = -(is + 4)H 2 + 3(2lu + is + 1)HF+[(is + 8)lu-6]F 2 

(JjUj — 

+ (u + 2)cm 2 + 3fc( ^ ± 4) . (2.15) 
2a z 

Equation (2.10) serves as initial condition for this dynamical system. 



3 Symmetry group analysis of basic 
dynamical equations 

We will use classical Lie group analysis [7] for symmetries of the system 
(2.14), (2.15) as a tool for finding exact solutions of these equations. Since 
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the equations are of first order, it is impossible to find all Lie groups of 
point symmetries that are admitted by this system (a symmetry condition 
cannot be split in the derivatives of unknowns to generate an overdetermined 
system, that can be solved and yield all the point symmetries, because these 
derivatives are expressed from (2.14) and (2.15)). Therefore, we make an 
ansatz for the form of symmetries that could be admitted by these equations 
under certain conditions. We note that equations (2.14) and (2.15)) are very 
close to those ones that admit scaling (or dilatational) symmetries and an 
obstacle to this is the mass term and/or curvature term (the one with k) 
in each equation. Indeed, a scaling symmetry group of transformations of 
independent and dependent variables is defined by 

a = X a a, F = X I3 F, H = X~<H. (3.1) 

It is generated by the infinitesimal generator 

x = aa i +pF w + " H m- (3 ' 2) 

Under transformations (3.1), equations (2.14) and (2.15)) become 

(2lu - 3)aH^- = {vu + 6)H 2 + A /3 " 7 2(z/ + A)uHF 

+ A 20J- 7) 2^ [{y + 6)uj _ 3]F 2 + A -2 7 2^ {vuj + 3)m 2 (3 3) 

+ A - 2 ^)A [^ + 2V + 3], 

(X 

-(2a; - 3)aH— = A 7_/3 - {v + A)H 2 + 3(2cj + v + 1)HF 
da 2 

+ X^'[{u + 8)cu - 6]F 2 + X- { ^'\u + 2)ujm 2 

+ A -( 2Q+ , +7 ) 3^ + 4) _ (34) 
2a 2 

The condition of invariance of equations (2.14) and (2.15)) under the Lie 
group of transformations (3.1) is that the transformed equations (3.3) and 
(3.4) should coincide with equations (2.14) and (2.15)), that is, all the A- 
dependence should vanish. 

The first obvious condition is 

7 = P, (3.5) 
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because it alone eliminates A in three terms, so that A will be present only 
in the terms with k and m. We have to consider several cases of constraints 
on k and m for vanishing the remaining dependence on A. 

1. Case 1 (generic). 

a + (3 = and a = 0, so that a = (3 = 7 = and hence the symmetry 
generator X = in (3. 2). There are no symmetries in the generic case. 

2. Case 2. 

k — 0, arbitrary m. Then (3 = 7 = and a is arbitrary, so it can be 
set to 1. The symmetry (3.2)becomes 

X = a Ta (36) 

3. Case 3. 

m = 0, arbitrary k. Then a + (3 = and, since 7 = (3, we can set 
/? = 7 = 1 and a = — 1 in (3.2), so that the symmetry becomes 

d d d . 

X = H 3H + F 3F- a Ta- < 37) 

4. Case 4. 

k = m = 0. Then /3 = 7 is the only condition and we have two 
inequivalent choices: 

(a) (3 = 7 = 0, a = 1 in (3.2) with the resulting symmetry 

X,=«|, (3.8) 

(b) /3 = 7 = 1, a = and the symmetry is 

To summarize, we have only one symmetry in cases 2 and 3 and two symme- 
tries in case 4. One symmetry is not enough for integrating a system of two 
first order equations, so in cases 2 and 3 we can only find invariant solutions 
[7], particular solutions of our equations. In case 4 with two symmetries 
we can find all solutions of our equations, that is, integrate our system in 
quadratures. 
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4 Case of flat spacelike sections: k = 



In this case we have only one symmetry (3.6) with the basis of invariants 
{H,F}. The equation expressed only in terms of invariants with solutions 
of the form F = F(H) can be obtained by dividing one of equations (2.14), 
(2.15) over the other one, with k — 0, with the result of the form dF/dH = 
G(H, F) which is a first order equation with no more known symmetries for 
m^O and hence it cannot be integrated. Though we cannot find the general 
solution in case 2, we still can search for invariant solutions with respect to 
the symmetry (3.6), which satisfy the conditions 

and produce a "static" solution {H,F} = constant. These H and F are 
roots of the algebraic equations 

{yu + 6)H 2 + 2(i/ + A)uuHF + + 6)u - 3}F 2 

+ ^(uuj + 3)m 2 = 0, (4.2) 

- (v + A)H 2 + 3(2c; + v + l)HF + [(u + 8)uj - 6]F 2 + (u + 2)um 2 = 0. 

There are four roots of this system of two quadratic equations. Two of them 
are real: 

H - ± 2m v^ F — ± =± F = - H 

^-uj(v 2 + Qv) - 12 ' yZ-Luiv 2 + Qv) - 12 " ^ 2 ' 

(4-3) 

where the expression under the square root in the denominators is positive 
for —6 < v-i < v < v\ < with the roots of the denominator = 

3 (-1 ± ^1-4/(3^)) and to > 10 4 . The latter range of values of v includes 

the interesting values Vdust = —3 and irradiation = —4. For these v and uj, two 
other roots of (4.2) are imaginary and will not be considered 

H = ± ^ ~ , F= T / , (4.4) 

- 17w + 12 V Quo 2 - 17u + 12 

where the roots of the denominator are uj\ = 3/2 and 0J2 = 4/3, so that 
6u 2 - 17u + 12 > 0. 
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The equations that define H and F 



a <p v 

— = H = constant, — = F = — H 
a '02 

are integrated to give time dependence of the solution 

a = a e Ht , = o ef m . (4.5) 

The solution (4.5) satisfies the initial condition (2.10) if arbitrary constants 
of integration O and ao are related by the formula 



*, = ±£J° J ( 7 6 t + T . (4-6) 
m V 3 V v\uj — 1) + 1 

where C is the coefficient of power law for the density (2.7). 



5 Massless case: rn = 

In this case we also have only one symmetry (3.7) with the basis of invariants 
{(p = aF,ip = aH} (indeed, Xip = 0, Xip = 0). Substituting F = (p/a and 
H = ip/a into (2.14) and (2.15) with m = 0, we obtain equations with the 
new unknowns if and ip: 

(2u - 3)a^^ = [(y + 2)cu + 3] (^ 2 + k) + 2{v + 



da 



+ ^[(z/ + 6V-% 2 , (5.1) 

-(2iu-3)a^ = ^(z/ + 4)(^ 2 + A;) + [4cu + 3(z/ + 2)]^ 

+ [(u + 8)oo-6] V 2 . (5.2) 

We can obtain only one equation expressed solely in terms of invariants by 
dividing one of these equations over another with the result of the form 
dip /dip = G((p,tp), but this equation has no more known symmetries for 
k 7^ and hence it cannot be integrated. Though we cannot find the general 
solution in case 3, again we can search for invariant solutions with respect to 
the symmetry (3.7) now, which satisfy the conditions 



9 



so that now ip and ip are constants independent of a. With these conditions, 
equations (5.1) form an algebraic system: 

\{u + 2)uo + 3](^ 2 + k) + 2(i/ + 4)u(p4> + + 6)u - 3]y? 2 = 0, 

(5.4) 

3(i/ + 4)(^ 2 + fc) + 2[4cj + 3(i/ + 2)]y^ + 2[(i/ + 8)cj - 6]y? 2 = 0, 
which again has four roots. Two real roots are 



(u + 2)VSk , VQk 



\/-2[(z/ 2 + 8z/)cj + 12cj + 6] ' ' v/-2[(z/ 2 + 8z/)cu + 12a; + 6] ' 

(5.5) 

that also yields F = ip/a and H = ip/a. Formulas (5.5) imply 

^t±V^ F = t±V H . (5.6) 

The expression under the square root in the denominators is positive for 
the values of v satisfying —6 < v<i < v < v\ < —2 where the roots of 
the denominator z/ 12 = —4 ± a/4 — Q/lo and to > 10 4 . Physically interesting 
values Vdust = _ 3 and v ra diation = ~ 4 again lie in this range, so the roots (5.5) 
are indeed real. Two other roots are (/2 = O=^F = O=^0 = O = constant, 
ip = ±\/—k. Thus, if k = 0, H = F = 0; if k = 1, H is imaginary; if k = — 1, 
ip = ±1, H — ±l/a and a = a ± t, that is a trivial solution. 

For the first two real roots we have a = ip = constant and so a = ipt + a . 
Then the relation <j>/<f> = ip/a yields time dependence of our solution 

= o (^ + a o )^ = M^t + ao) {u+2)/2 , a = ^jt + a , (5.7) 

where we have used (5.6). The solution (5.7) satisfies the initial condition 
(2.10) if 0o is expressed in terms of other constants as follows 

2v/OJ 

00 = ±- 



2 



a/3(^ 2 + fc) + 6^ + 2ujif 

= ± , ^ (5.8) 

y/{u + 2)[(z/ + 2){2u + 3) + 12]^ 2 + 12k 

where C is again the coefficient of power law for the density (2.7). 
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6 Massless case with flat spacelike sections: 
k = m = 



In this case we have two commuting symmetries (3.8) and (3.9) which is 
enough for the complete integration of our equations in quadratures that 
will yield their general solution. We apply first the symmetry generator 
Xi = ad /da with the invariants H and F. A solution expressed solely in 
terms of invariants has the form F = F(H). Dividing equation (2.15) over 
(2.14) with k = m = 0, we obtain a single first order equation for F(H): 

dF [(z/ + 8V-6]F 2 + 3(2cu + z/+l)#F+(3/2)(z/ + 4)ff 2 

(6.1) 



dH (2w/3) [(u + 6 V - 3]F 2 + 2{y + A)ujHF + {voo + 6)H 2 ' 

This equation still admits one symmetry X 2 = Fjfp + with the in- 

variant G = F/H (and also a, not contained explicitly in (6.1)), so that a 
new independent variable is H and a new unknown is G = G(H). In the new 
variables {H,G} equation (6.1)) becomes the one with separated variables 
H and G 

dH 2co[(u + Q)lo - 3}G 2 + Q(u + 4)wG + 3(uco + 6) 



if {2[(z/ + 6)cj - 3]G + 3(i/ + 4)}(2cjG 2 + 6G + 3) 

where the cubic polynomial in G in the denominator was factorized. Inte- 
grating both sides of this equation in if and G and getting rid of logarithms, 
we obtain explicitly the function H(G) 

H = H ({2[(u + 6)u - 3]G + 3(i/ + 4)} i/ (2cjG 2 + 6G + 3) 3 )~ 6 x (6.3) 



2luG 2 + QG + 3 




z/ + 3)/ tan 



2[(z/ + 6V-3]G + 3(z/ + 4) 
where we have introduced the following constants 



2uG + 3 
v/3(2cj - 3) 



(z/ + 6V-3 2^3(2^-3) 

6 



(z/ + 6) 2 cj-6(z/ + 5) ' (z/ + 6) 2 cj-6(z/ + 5) ' 

= 3(1/ + 4) 

9 (z/ + 6) 2 cu-6(z/ + 5) ' 1 j 

Using (6.3), from the definition of G we also have F as a function of G: 
F = GH(G). 
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The dependence G(a) can be determined from (2.14) with F = GH(G) 
at k = and m = 0: 



da 
a 



3(2cu-3) 



dH 



x 



2u[(v + 6)cj - 3]G 2 + 6(1/ + 4)wG + 3(z/cj + 6) " // 
where dH/H is expressed from (6.2). The resulting equation is 



da 
a 



-6(2u - 3)dG 



{2[(v + 6)u - 3]G + 3(u + 4)}(2cuG 2 + QG + 3) 



(6.5) 



Integrating both sides of (6.5) and eliminating logarithms, we obtain an 
explicit dependence a(G) 



a 



2uG 2 + 6G + 3 



{2[(z/ + 6V-3]G + 3(z/ + 4)} s 



x exp < / tan 



-i 



2uG + 3 
v/3(2cj - 3) 



(6.6) 

Formulas (6.3) and (6.6) together with F = GH(G) yield a parametric rep- 
resentation, with a parameter G, of the general solution H(a), F(a) to equa- 
tions (2.14) and (2.15) with k = m = 0. 
From 616 — F — GH we have 



d6 
J 



GHdt = G 



-6(2u - 3)GdG 



da 

~~a ~ {2[(i/ + 6)w - 3]G + 3(i/ + A)}(2ujG 2 + 6G + 3) 



where we have expressed da/a from (6.5). Integrating both sides of this 
equation and eliminating logarithms, we obtain the explicit dependence 6(G) 



x exp 



2[(z/ + 6V-3]G + 3(z/ + 4) 
V2ujG 2 + 6GT3 

2luG + 3 



v/3(2cj - 3) 



(6.7) 



and because of the known dependence a(G) from (6.6) we have determined 
6(a) in a parametric representation. In order to satisfy the initial condition 
(2.10), the constants in the formulas (6.3), (6.6), and (6.7) should be related 
as follows 

u/2 



(6.8) 
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To determine time dependence of the solution, we use a/a = H to obtain 

/da 

and then express da/a from (6.5) and H from (6.3). The integrand simplifies 
if we introduce the new parameter g = (2ujG + 3)/ ^3{2u — 3) with the 
resulting integral 



t-t n = 



3(2cu-3) (2 



Ho 



{[(u + Q)co - 3}g - ^3(2uj - 3)}2("+3)6-i 



(g2 + l)(v+3)b 



X 



exp [— {y + 3)/ tan 1 g\dg, 



(6.9) 



that gives us explicitly functions t(g) and hence t(G) and implicitly the time 
dependence G(t). Then from the functions a(G) and 4>(G), determined by 
(6.6) and (6.7) respectively, we know the time dependencies a(t) and <j)(t) in 
a parametric representation with the parameter g in terms of the quadrature 
(6.9). 

For dust, v = u dust = —3 and the integral (6.9) simplifies substantially 
with the result 



t-U, 



{2uo - 3) 



In 



3(w-l)g- ^3(2lu - 3) 



(6.10) 



(to - I) H 

that can be inverted easily to yield g(t) and so an explicit expression for G(t): 
1 



G = 



2u(u - 1) 



2cj-3 



exp 



(u; 



{2uo - 3) 



H (t - 1 ) 



00 



(6.11) 



For radiation, v = ^radiation = ~4 and the simplification is not enough for 
evaluating the integral in (6.9). 



7 Conclusion 

We have shown that, in the model of cosmology we are considering, the 
equation of state does not affect the expansion law of the universe. Thus, for 
the m = case the expansion of the universe, following the invariant solution 
(5.7), is always linear. What is interesting is that according to (5.5) this 
solution is nontrivial only for k=l i.e. a closed universe. It is also interesting 
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to note that the equation of state which determines the coefficient v affects 
the rate of expansion. This case may be relevant for the early universe where 
the curvature term cannot be neglected. (5.7) shows that with ao = 0, for 

t « mr l 

0Z/ + 21 z/ + 2 

4, = — 1»— m ' <"» 

so that m\<f>\ << \<p\ and for early times the mass term can indeed be ne- 
glected. After the universe expands to a large size, the curvature term can be 
neglected and in this case the k=0 solution becomes important. (4.3), (4.5) 
show that the universe expands exponentially for all equations of state. The 
equation of state is important just for the value of the Hubble constant. This 
behavior we identify with the dark energy or cosmological constant which is 
observed today. 
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